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RESEARCH MEMORANDUt-I 


STATIC AEROELASTIC PHENOMENA 
OF M-, V-, AND A-WINGS 

By Franklin W. Diederich and Kenneth A. Foss 


SUMMARY 


Spanwise lift distributions, lift coefficients, spanwise centers of 
pressure, shifts in aerodynamic center, coefficients of damping in roll, 
aileron rolling-moment coefficients, and rates of steady roll per unit 
aileron deflection have been calculated for nine M-, W-, and A-wings, 
as well as for comparable ordinaiy sweptforward, unswept, and sweptback 
wings. Although the calculations are too specific to permit any quanti- 
tative conclusions which are generally applicable, certain qualitative 
conclusions are drawn conceraing the plan forms most suitable from the 
aeroelastic point of view. In general, there is reason to believe that 
certain M and W plan forms exist which are superior aeroelastically and 
structxirally to ordinary swept wings. 


INTRODUCTION 


The use of M- and W-wings has been suggested as a means of allevi- 
ating the static aeroelastic problems of swept wings, such as the shift 
of aerodynamic center and the loss of lateral control. The obvious, but 
not necessarily most economical, remedy for these static aeroelastic 
difficulties consists in stiffening these wings. The advantage of an 
M- or W-wing over an ordinary swept wing is that, inasmuch as in the M- 
or W-wlng the over-all effects of bending and torsion deformations tend 
to oppose each other, a more flexible structure may be acceptable. The 
concept of alleviating static aeroelastic difficulties by such a 
balancing process is not new. In reference 1, for instance, means are 
discussed for achieving this balance organically with a swept wing, and 
in reference 2 an artificial balancing device (an auxiliary lifting sur- 
face mounted on a boom of the tip of a sweptback wing) is analyzed. 

In the appendix of the present paper a method based on those of 
references 3 and 4 is presented for analyzing the static aeroelastic 
phenomena of M-, W-, and A-wings, such as change in aerodynamic loading, 
lift-curve slope, and aerodynamic- center shift due to aeroelastic effects, 
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divergence j loss of lateral control, and change in damping in roll due 
to aeroelastic effects, as well as aileron reversal. 

This method has "been used to calculate the foregoing aeroelastic 
phenomena for three M-, three W-, and three A-wings (including one 
inverted A-wing) . The results of these calculations are discussed with 
particular regard to the problem of selecting optimum plan foimis for 
the minimization of the adverse effects of static aeroelastic phenomena. 
Similar calculations have also been made for a sweptbach, a sweptforward, 
and an unswept wing to afford a basis of comparison of the static aero- 
elastic phenomena of the M-, W-, and A-wings with those of the more con- 
ventional wings. The calculated aeroelastic phenomena of the M-, W-, 
and A-wings are discussed in the light of these comparisons. 


SYMBOLS 


A aspect ratio, b^/s 

a distance of section aerodynamic center from leading edge, 

fraction of chord 


a 


Aa 


position of wing aerodynamic center measured from leading 


edge of mean aerodynamic chord 


A 

U ‘ ct y 


wing aerodynamic-center shift 



b wing span 

bg span of both ailerons 

b ' span of exposed wing (b - w) 

Cg bending moment (semispan rolling-moment) coefficient, 

l4M^yfcSb 

lift coefficient, L,/qS 




a 


lift curve slope 


Cj rolling-moiTient coefficient for a linear antisymmetrical 

angle-of-attack' distribution with a tip angle of one 

radian, -C 2 ^ 

P 
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CZp 


c 

•-a 
= 2 


c 


c 


d 

El 


e 


®1 

®2 

GJ 

L 

2 

M 


damping- in- roll coeffic-ient 

rolling-moment coefficient due to a unit aileron 
deflection 

pitching-moment coefficientj Mp^gSc 

chord of wing (measured parallel to air stream) 

chord of aileron 

section lift coefficient, 2/ qc 

chord at airplane center line 


average chord, s/b 


mean aerodynamic chord. 



span of part of wing, fraction of b/2 

bending stiffness in planes perpendicular to elastic 
axis 


distance of elastic axis from leading edge, fraction of 
chord 

dimensionless moment arm (e - a) 


distance of center of pressure due to aileron deflection 
behind elastic axis, fraction of chord 

torsional stiffness in planes perpendicular to elastic 
axis 

lift on total wing span 

lift per 'unit distance along span 

accumulated bending moment about axes parallel to air 
stream 
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accumulated, bending moment about axes perpendicular to 
elastic axis 

accumulated single-wing rolling moment about fuselage 
center line 

total wing pitching moment about quarter-chord point of 
mean aerodynamic chord, positive nose up 


wing-tip helix angle due to roll 


root-rotation constants defined in equations (33) # (3^)^ 
and ( 36 ) , respectively 


dynamic pressure 


dimensionless dynamic-pressure parameter, 

(GJ)r 


wing area 

accumulated torque about axes perpendicular to plane of 
symmetry 

accumulated torque about elastic axis 
running torque in planes parallel to air stream ( 
width of fuselage at wing root 
distance defined in figure 1 

streamwise distance of a section aerodynamic center aft 
of an unswept reference line through the quarter- chord 
point of the mean aerodynamic chord 

lateral ordinate (see fig. 1) 

dimensionless lateral ordinate, 

b/2 

dimensionless position of lateral center of pressure 



a 


°6 


r 

5 

A 

\ 


angle of attack measured in planes parallel to air 
stream 

angle of attack equivalent to unit aileron deflec- 
dc 7 /d6 

tion, 

dc2/da 

local dihedral angle (in a plane through elastic axis) 
due to wing deformation along elastic axis 

aileron deflection measured in planes parallel to air 
stream 

angle of sweephack at elastic axis 

taper ratio, c^jc£ 


(p angle of twist in planes perpendicular to elastic axis 

Subscripts: 

B at point of wing break (point of spanwise discontinuity 

of angle of sweep) 

£. at center line of airplane 


D 

f 

g 

i 

o 

R 

r 

s 

w 


at divergence 

that portion of wing covered by fuselage 

geometric (built in or due to airplane attitude) 

inner part of wing, not including that part covered by 
fuselage 

outer part of wing, from wing break to wing tip 

at aileron reversal 

wing root (located at Intersection of elastic axis and 
fuselage side) 

structural (due to structiaral deformation) 

wing alone (not including that portion covered by 
fuselage) 

pertaining to rigid wing (q = O) 


0 
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Matrix Notation: 

Note: Specific matrices are defined where they first occur. 


II II 
[] 
L1 

{} 

LJ 


rectangular matrix 
square matrix 
diagonal matrix 
column matrix 
row matrix 


DESCRIHTIOK OF THE CALCULATIONS 
Method of the Calculations 


The method used to perform the calculations of static aeroelastic 
phenomena of M-, and A-wings is presented in the appendix to this 
paper j it is based on the method of references 3 and k and consists, 
like those methods, in integrating by means of numerical and matrix 
techniques the differential equations which describe the static aero- 
elastic phenomena. 


A.ssumptions 

The spanwise lift distribution is assumed to be given by suitable 
aerodynamic influence coefficients and the local centers of pressure 
of the lift due to angle of attack and due to aileron deflection are 
assumed to be invariant with angle of attack and aileron deflection. 
Both of these assumptions imply small angles’ of attack and aileron 
deflection. 

A straight elastic axis is assumed to exist in both parts of the 
wing, and the wing is assumed to be mo\mted flexibly at an effective 
root perpendicular to the elastic axis through the Intersection of the 
elastic axis and the fuselage (see fig. l) so that; the root triangle 
imparts rigid-body rotations to the wing, the rotations being propor- 
tional to the root bending moment and the root torque. On the other 
hand, the outer part of the wing is assumed to be attached rigidly to 
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the inner part, so that, if the inner part were rigid, no rigid-hody 
rotations would he imparted to the outer part of the wing. All defor- 
mations beyond those due to the rigid-hody rotations imparted hy the 
root triangle are then assumed to he given hy the elementary theories 
of bending and of torsion along the elastic axis. 

The angle between the aileron and the wing is assumed to he 
constant along the span of the aileron. This assumption implies that 
the aileron and wing twist the same amoimt. 


Scope of the Calculations 

The M-, W-, and A-wing plan forms for which calculations have been 
made are listed in table 1 as wings 1 to 9. Wings 1 to 3 are M-wingsj 
wings 4 to 6 are W -wings; wings 7 and 8 are A-wings and wing 9 is an 
inverted A-wing. For the sake of comparison, calculations have been 
made also for three conventional plan forms - a sweptforward wing, an 
unswept wing, and a sweptback wing - listed in table 1 as wings 10, 11, 
and 12, respectively. All wings have a taper ratio of 0.5j all have 
angles of sweep of either zero or ±45°, and all have an aspect ratio 
of 6. Three values of the spanwise position of discontinuity in sweep, 
hereinafter referred to as the "break,” are included in this series of 
plan forms, namely, y* = 0.3, 0.5, and 0.?. All wings were considered 
to be mounted on a fuselage of width equal to 0.1 of the span. 

For all plan forms, symmetrical lift distributions were calculated 
for one subsonic and one supersonic flow condition at values of the aero-- 
elastic parameter q equal to 3*0, and for most plan forms for q = 6.0 
as well. Lift distributions were calculated for unit geometric angle of 
attack across the span, for linear antisymmetric geometric angle of 
attack with unit angle at the tip; and for unit effective angle of attack 
due to the deflection of an outboard aileron. These lift distributions 
were Integrated to obtain total lifts, rolling moments, and positions 
of the wing center of pressure. 

For subsonic speeds the lateral- control parameters were calculated 
for 20 -percent- chord, 50-percent- span outboard ailerons, with some addi- 
tional calculations for all-movable wing tips extending over the outer 
30 percent of the semispan (which may be considered to be 100 percent- 
chord ailerons) . For supersonic speeds lateral- control calculations 
were made for 20-percent-chord, 30 -percent- span outboard ailerons, with 
some additional calculations for 50-percent-span ailerons. (See table 2.) 


Basic Data 

The spanwise stiffness distributions used in this paper are given 
in figure 2. These stiffnesses are based on the constant-stress analysis 
in reference 1 (with modifications occasioned by the wing break) , except 
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that the stiffnesses were taken proportional to the fourth power of the 
chord from y = 0.7 to 1.0 so that there would he finite stiffness at 
the wing tip. Other structural parameters (including the root-rotation 
constant) are given in table 1; the values of the stiffness ratio given 

A , , > 

( = 0.79^1 and of the elastic-axis location (e = 0.45) are typical 

\(Sl)r / 

of wings of conventional thin- skin construction, having a thickness 
ratio of about 10 to 12 percent. 

For subsonic speeds aerodynamic influence coefficients were calcu- 
lated by the method of reference 5; the manner in which this method was 
modified to apply to M-, W-, and A-wings is discussed in the appendix. 
The rigid-wing subsonic lift distributions used to compute the aero- 
dynamic influence coefficients were taken from reference 6, in which 
they were calculated for incompressible flow by simplified lifting- 
surface theory. 

For supersonic speeds strip theory was used because no suitable 
means (such as aerodynamic influence coefficients) were available for 
calculating lift distributions for angle-of-attack distributions which 
are not Initially known, although lift distributions can be calculated 
for any given angle-cf-attack distribution by linearized supersonic 
theory. (See, for Instance, the methods of refs. 7 and 8.) The devel- 
opment of such coefficients solely for the purpose at hand was not con- 
sidered justified principally because M-, W-, and A-wings are intended 
primarily for flight at subsonic and transonic speeds. Also, the 
results of the calculations of the present paper can be interpreted as 
flexibility corrections to the aerodynamic characteristics of the rigid 
wings; if the corrections are relatively small the corrected results 
will be relatively insensitive to the assumptions made in calcialating 
the corrections. The resulting flexible-wing characteristics are then, 
of course, no better than rigid-wing characteristics to which the cor- 

I rectlons are applied. Inasmuch as in the present paper the aeroelastic 

• increments are of primary Interest, the rigid-wing characteristics were 

• estimated by strip theory for the sake of simplicity. 

The subsonic local aerodynamic-center positions were also taken 
from reference 6; and the corresponding dimensionless section moment 
arms, e^^ = e - a, are plotted in figure 3- For convenience, the local 

aerodynamic centers for supersonic flow were assumed to lie along the 
45-percent-chord line, so that they coincide with the elastic axis; 
thus e]_ is 0 along the entire span. 

For subsonic speeds the values of the dimensionless section moment 
arm due to aileron deflection, eg, were calculated from the two- 

dimensional values by assuming that the difference between the two and 
three-dimensional lifts acts at the section aerodynamic center as 


NACA EM L52J21 


described in references 4 and 9- The values of &2 obtained in this 
manner are plotted in figure 3- For the wings with rotating tip panels 
the values of eg were assumed to be equal to ( -e]_) . For supersonic 
speeds the center of pressure due to aileron deflection was assumed to 
be at 90 percent of the chord; consequently, a value of eg = 0.4-5 was 
used for all wings. 

The quarter-chord point of the mean aerodynamic chord is used as a 
reference for pitching moments. Unlike the case of an ordinary wing the 
longitudinal location of the mean aerodynamic chord does not coincide 
with that of the chord at the station which corresponds to the centroid 
of area of the wing. An expression for the distance of the quarter- 
chord point of the mean aerodynamic chord rearward of the intersection 
of the quarter-chord line and the plane of symmetry is given in equa- 
tion (lO) of reference 6; however, attention is called to the fact that 
in reference 6 the angle of sweepback refers to the quarter-chord line 
rather than the elastic axis. 


Results of the Calculations 


Spanwise lift distributions .- The spanwise lift distributions for 
the nine M-, W-, A-, and three ordinaiy wings are presented in figures 4 
to 15 for two or three dynamic pressures including 0 (the rigid-wing 
case) and for subsonic as well as for supersonic speeds. The top parts 
of the figures show the lift distributions due to a unit airplane or 

cci 


root angle of attack represented by the coefficient 

cC 


the lift 




coefficient being that of the given wing for q = 0. The bottom 

parts of the figures show the lift distributions due to unit effective 

cci 


aileron deflection (ccg6) represented by the coefficients 


that 


cC- 




c 


is, the loading coefficients 
of the rigid wing. The coefficient 


divided by the damping coefficieht 
cci 


cC 


Ido 


the -product of the coefficients 




and 


can also be construed as 


the first of these 



coefficients represents the loading coefficient per unit rolling moment 

/ cc^ \ 

of the rigid wing and is analogous to the coefficient [— 1 used for 


cC 
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the lift dlstrlDution due to angle of attackj whereas the second coef- 
ficient is the wing-tip helix angle per unit effective aileron deflec- 
tion of the rigid wing. 

The calculations nade in reference 6 which form the basis for the 
aerodynamic information used for the subsonic calculations in the pres- 
ent paper pertain to a Mach number of 0, and hence, so do the calculated 
results. However, the lift distributions may be expected to be substan- 
tially unchanged (for small angles of attack) throughout the subsonic 
region and, except locally near the fuselage and the wing break and 
except for the unswept and A-wings through the transonic region as well. 
The lift distributions for supersonic speeds were estimated by strip 
theory and are, therefore, independent of Mach number. 


If the rigid-wing lift distributions are independent of Mach number 
so are the flexible-wing lift distributions. However, inasmuch as the 
rigid-wing lift-curve slope ^Loq enters into the definition of q the 

lift distributions for a given value of q correspond to different 


values of 


if C 




°0 


changes with Mach number. 


Aerodynamic parameters associated -with the spanwise lift distri- 
bution .- The quantities Cl/Clq> , a, 

and pb/2V are presented in table 2 and the quantities a, 

C2g^C7g , and pb/2V are also plotted in figures l6 to 22 against, the 

dimensionless dynamic pressure q for several of the M-, V7-, and A-wings 
as well as for the unswept and sweptback wings. The values of 

and a given in table 2 were extrapolated to the large values of q 
represented in figures l6 to 22 by the use of the approximate formulas 
in reference 1. The values of the coefficients C^, and C 2 g for 

dynamic pressures other than 0 can be determined from the ratios 
C 2 ^^C 2 p , and > since the rigid-wing values ^^p ' 


=0 


are presumably known. The values of Cjj and C 2 are given 


°0 




in table Ij CLq can be obtained from value of the air- 

plane angle of attack, and 02 ^ is equal to the product of (pb/2V)Q 
and C 2 j • (The values of (pb/2V)Q are the values of pb/2V given in 


table 2 for 


= 0 .) 
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l?he rigid-wing values of and given in table 1 for sub- 

sonic speeds are those calculated in reference 6 for M = Oj in principle 
they can be corrected for subsonic compressibility effects by the three- 
dimensional Prandtl-Glauert rule, but calculations must be available for 
M-, W-, and A-wings with many different sweep angles and aspect ratios 
before this collection can be effected. The correction for can 


be used for well; within this approximation 

unaffected by compressibility at subsonic speeds. 


- 


-"0 
s then 


The rigid-wing values of and given in table 1 for super- 

sonic speeds are estimated on the basis of the Acheret theory for M = 2 
and are intended for qualitative comparisons only. For quantitative 
purposes they can be calculated by linearized supersonic theory not only 
for M = 2- but for any supersonic Mach number which is not too large 
nor too close to 1. For ordinary wings the results of such calculations 
are presented in references 10, 11, 12, and 13^ for instance. 

In considering figures l6 through 22 the fact should be kept in 
mind that the abscissa is subject to compressibility effects to the 
extent that, as previously mentioned, for a given value of q a change 
in *^100 implies a change in q. Apart from this effect the resiolts 

presented in these figues are independent of compressibility effects 
provided the lift distributions (within a given speed region) are sub- 
stantially unaffected by Mach number. 


The dynamic pressures at divergence and at reversal .- The values 

of the dimensionless dynamic pressure q at divergence and at aileron 
reversal are given in table 2. From these values the corresponding 
values of q can be calculated from the definition of q. In cases 
for which the lowest (in absolute magnitude) dynamic pressures required 
to diverge the wing were found to be negative, the next higher critical 
dynamic pressures were calculated by using the method outlined in this 
paper; these values are also presented in table 2. For the wings with 
rotating tips the lowest (in absolute magnitude) dynamic pressures 
required to reverse lateral control were foimd to be complex. 


DISCUSSION 

Comparison of the Aeroelastic Properties of the Various Wings 


the 


Spanwise lift distributions .- The rigid-wing lift distribution of 
unswept wing is approximately elliptical at subsonic speeds, and the 
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effects of sweepforward and sweepback on the spanwise lift distribution 
are to shift the load inboard and outboard, respectively. As pointed 
out in reference 6, the rigid-wing lift distributions calculated therein 
for these compoionded plan forms at subsonic speeds are similar to those 
that could have been estimated qualitatively from the knowledge of the 
characteristics of the lift distributions of the ordinary swept wings 
of which the compounded plan forms may be considered to be composed': In 

the case of the M-wing represented in figure h, for instance, the inner 
part of the wing behaves aerodynamically like a sweptforward wing with 
the characteristic peak in the lift distribution near the plane of sym- 
metry, (see fig. 13 ), whereas the outer part behaves like a sweptback 
wing with the characteristic loading up of the wing tip (see fig. 15)- 

As previously mentioned, for supersonic speeds strip theory was 
used to estimate the rigid-wing lift distributions as well as the aero- 
elastic increments to these distributions. The rigid-wing lift distri- 
butions are therefore identical for al.1 wings. 

At subsonic speeds the effect of aeroelasticlty on the xinswept wing 
is to increase the lift, particularly near the tip, in the symmetric 
case, and to decrease the lift in the aileron-deflected case. (See 
fig. ll<-.) At supersonic speeds aeroelasticlty has no effect on the 
symmetric lift distribution because the center of pressure was assumed 
to be in a position which coincides with the elastic axis. The decrease 
in lift in the antisymmetric case is quite pronounced, however, due to 
the fact that the moment arms B 2 (or, more to the point, the sums of 

the moment arms eg + are relatively large . 

The symmetric lift distributions on the sweptforward wing (fig. 13) 
at subsonic speeds exhibit an even larger increase in lift due to aero- 
elastic effects than do those of the unswept wing. At supersonic speeds 
there is also a large increase in lift on the sweptforward as compared 
to the unswept wing; this increase is due entirely to the bending of the 
wing, inasmuch as the moment arm e^ is zero as for the unswept wing. 

The lift due to aileron deflection is Increased as a result of aero- 
elasticity because in the case of the sweptforward wing the bending 
effects which tend to increase the lift due to aileron deflection pre- 
dominate over the torsion effects which, as in the case of the unswept 
wing, tend to decrease this lift. 

On the sweptback wing (fig. 15) the effect of the bending defor- 
mations also predominates over the torsion deformations but causes a 
decrease in lift in the symmetric case and augments the effect of the 
torsional deformations in the aileron- deflected case as q increases 
to produce a large loss of lift. 

The effects of aeroelastic action on the lift distribution of the 
compounded plan forms are qualitatively as may be expected from a 
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knowledge of the aferoelastic effects on the lift distritrations on the 
constituent parts of the wing. In the case of the M-wing 1 (fig. 4), 
for instance, the large sweptback outer part of the wing results in 
aeroelastic effects which are similar to those of a sweptback wing in 
that they decrease the lift both in the symmetrical case and in the 
aileron-deflected case, both at subsonic and supersonic speeds. The 
aeroelastic effects on the lift distribution on the sweptforward inner 
part of the wing are in the opposite direction but not large enou^ to 
result in an increase in lift but merely to decrease the loss in lift 
caused on the inner part of the wing by the aeroelastic action of the 
outer part of the wing. 

The aeroelastic effects on the lift distribution of wing 2 ( see 
fig. 5) are similar to those on the lift distribution of wing 1, but 
due to the large relative size of the sweptforward inner part of the 
wing an increase in lift is actually noted in the symmetric case on the 
inner part of the wing. As a result of this Increase there is a tend- 
ency for the aeroelastic effects on the symmetric lift distribution to 
cancel. If the inner part of the wing were sli^tly larger still, the 
lift and center of pressure would probably be substantially unaffected 
by aeroelastic action. The loss in lateral control due to aeroelastic 
action is less than that of wing 1, but still quite large. 

In the case of wing 3 (fig* 6) the position of constant lift and 
center of pressure has been passedj the aeroelastic characteristics of 
the large sweptforward inner part of the wing dominate the aeroelastic 
behavior of the wing in the symmetric case, although in the supersonic 
case the opposite aeroelastic characteristics of the sweptback outer 
part are sufficiently large to cancel the increase in lift resulting 
from the aeroelastic behavior imparted to the whole wing by its inner 
part, at least near the tip. In the aileron-deflected case the aero- 
elastic behavior of the inner part dominates that of the outer part and 
results in an increase in lift at subsonic speeds. At supersonic speeds, 
however, the outer part of the wing dominates the inner part to the 
extent that the loss in lift is only slightly smaller than that noted 
for wings 1 and 2. 

The W-wings k and 5 (see figs. 7 and 8) have large sweptforward 
outward parts which completely dictate the aeroelastic behavior of the 
entire wings; the sweptback inner parts are capable only of reducing 
slightly and locally the Increase in lift Imposed everywhere on the wing 
as a result of the aeroelastic action of the outer part of the wing. 

The W-wlng 6 (see fig. 9) is close to an over-all aeroisoclinlc 
condition, that is, a condition of over-all cancellation of the effects 
of bending and torsion deformations. In the 'subsonic symmetric case the 
lift, and in the supersonic symmetric case the center of pressure, are 
substantially unaffected by aeroelastic action as a result of the balance 
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between the aeroelastic tendencies of the inner and outer parts of the 
wing. There is an outboard shift of the center of pressure in the sub- ^ 

sonic case, as for wings A and 5> a decrease in lift in the super- 
sonic case, which is opposite to the behavior of wings 4 and 5- In the y, 

antisymmetric case, both at subsonic and supersonic speeds, the aero- 
elastic behavior of the large sweptback inner part of the wing dominates 
that of the whole wing and results in a loss in the lift due to aileron 
deflection. 

In the case of the A-wing 7 (see fig. lO) as in the case of 
W-wings 4 and the behavior of the outer part dominates that of the 
sweptback inner part except in the supersonic symmetric case, for which 
the outer part is largely inactive as far as aeroelastic behavior is 
concerned so that the aeroelastic behavior of the sweptback inner part 
resulus in a small decrease in lift. On the other hand, in the case of 
A-wing 8 (see fig. ll) the large sweptback inner peirt of the wing domi- 
nates its aeroelastic behavior. The aeroelastic action of the unswept 
outer part only serves to reduce the resulting loss of lift locally to 
a small extent. Similarly, the large sweptforward inner part of the 
inverted A-wing 9 (see fig. 12) largely dominates the aeroelastic 
behavior of that wing, except that at supersonic speeds the twist of 
the outer part resulting from the large moment arm 69 of the lift due 
to aileron deflection is so large that it overshadows the bending effects 
of the inner part of the wing and resxilts in a small loss in the lift due 
to aileron deflection. 

t 

In general, the effects of aeroelasticity on the spanwise lift dis- 
tributions may be seen to be much less for certain compounded plan forms 
(wing 6, for example) than for ordinary swept wings. 

Inasmuch as strip theory was used in the calculations for super- 
sonic speeds the results for supersonic speeds cannot be expected to be 
as accurate as those for subsonic speeds. If more realistic values of 
the lift distributions are desired for supersonic speeds the increments 
due to aeroelastic action shown in figures 4 to 15 can be applied to 
rigid-wing lift distributions calculated by linearized supersonic theory. 

These increments are probably quite accurate, because the integrating 
matrices used in the present paper have the effect of rounding off any 
lift distribution to which they are applied; the aeroelastic effects 
would have been overestimated slightly if strip theory had been used 
rigorously. 

Aerodynamic parameters associated with the lift distributions .- The 
lift and aerodynamic center are determined by the symmetrical spanwise 
lift distributions; similarly, the rolling-moment coefficient due to 
aileron deflection and wing-tip helix angle per unit aileron deflection 
are determined by the corresponding antisymmetric lift distributions. 

The effects discussed in this section are therefore direct consequences 
of those discussed in the preceding section. 
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The effects of aeroelasticity on some of the aerod 3 mamic properties 
of the \inswept wing are shown in figure 21. They result in an increase 
in lift and a loss in the rolling power and rate of roll at subsonic 
speeds; the rate at which lateral control is lost, slowly at first and 
then more rapidly, is typical of wings for which positive, 

as discussed in reference 9- At supersonic speeds the lift is unaf- 
fected, hut the losses in lateral control and rolling velocity are even 
greater than at subsonic speeds. The rate at which, the control is lost 
with increasing q is constant, a phenomenon typical of wings with 
infinite q^. 

Although the lateral control of the sweptforward wing is improved 
by aeroelastic action, the increase in lift and the outboard shift in 
center of pressure are so large (see table 2) as to make this type of 
plan form undesirable. The sweptbaclc wing 12 (fig. 22), on the other 
hand, experiences a loss in lateral control which is even greater than 
that of the unswept wing. The rate at which control is lost is rapid 
at first, then slower; as is typical of wings with negative. 

The sweptback wing loses some lift and its center of pressure moves 
inboard as a result of aeroelastic action. This movement of the center 
of pressure is accompanied by a shift of the aerodynamic center forward 
and, hence, a loss in the static- stability margin. The loss of control 
and the shift of the aerodynamic center are disadvantages of the swept- 
back wing from the aeroelastic point of view and the aim of this anal- 
ysis is to determine whether there are compounded plan forms which are 
substantially superior to the sweptback wing in this respect. 

As noted in the preceding section, there are among the nine com- 
pounded plan forms considered some wnich exhibit little or no loss in 
lift due to aileron deflection and little shift in spanwise center of 
pressirre as a result of aeroelastic actions. (Of course, the fore-and- 
aft movement of the center of pressure varies with the spanwise shift 
in a more complicated manner than in the case of the ordinary swept 
wings as a result of the complicated geometry of the compounded wings.) 
As shown in figures l6, l8, and 20, the aerodynamic-center shift of 
wings 1, 6, and 8 is smaller than for that of the sweptback wing 
(fig. 22), and the shift for wings 2 and 7 (figs. 17 and 19) is practi- 
cally nil. However, the loss of lateral control of wings 1 and 2 is 
only slightly less, and that of wing 8 is actually slightly greater, 
than that of the sweptback wing. At low dynamic pressures (q less 
than 6) wings 6 and 7 siiffer relatively little loss of control. These 
five wings are typical of the others, as may be seen from table 2, 
except that the other wings actually experience a gain in lateral con- 
trol. However, this gain is purchased at the price of greater shift of 
aerodynamic center (all but wing 5) or low divergence speed (particular! 
wings 4, 5^ and 9) with the resulting tendency to general aeroelastic 
instability. Thus, although some of the compounded plan forms have 
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generally ncre favorable aeroelastlc characteristics than does the 
sweptback wing^ some show little improvement and others are actually 
inferior to the sweptback wing from the aeroelastic point of view. 

Dynamic pressures at divergence and at reversal .- For ordinary 
unswept wings the dynamic pressure is usually primarily of interest as 
a reference quantity which serves as an index for the severity of static 
aeroelastic pheonomena. Only for sweptforward wings does the divergence 
speed have any physical significance, and for these wings it is likely 
to be so low as to rule out the use of the wings because stiffening the 
wings would require a prohibitive amount of structural material. On 
the other Znand, the dynamic pressure required to diverge sweptback wings 
is negative, so that its only significance is as a reference parameter. 

The significance of the dynamic pressures listed in table 2 can be 
assessed by the fact that a sweptfo27ward wing is likely to diverge at 
relatively low dynamic pressures. (See refs. 3 1^, for instance.) 

Therefore, a value of q = 6 may be expected to be attained by a 
fighter-type airplane at about Mach number 1 at low altitudes. This 
value of q is seen to be close to the value of divergence for some 
of the W-wings and the inverted A-wing; the sweptforward wing would 
diverge at an even lower value of q. Similarly, at a value of q = 6 
several of the M-wings would experience aileron reversal at supersonic 
speeds. However, more definite statements cannot be made unless the 
physical parameters that enter into the definition of that is, the 
dynamic pressure of operation as well as the geometric and structural 
properties of a given wing, are known. 

An interesting use of the dynamic pressure as a reference param- 
eter was pointed out in reference 1. As the sweep of an ordinary wing 
is var,ied from unswept to sweptback, or as the stiffness-ratio GJ/EI 
or the elastic-sjcis location of a sweptback wing is varied, the dynamic 
pressure required for divergence goes to infinity and then reverses sign 
at a particular combination of structural, geometric, and aerodynamic 
parameters. For this combination of parameters the bending and torsion 
deformations lead to forces which tend to cancel eacii other; in other 
words, the aerodynamic loads give rise to defoimiations which do not give 
rise to any further aerodynamic loads. This phenomenon is referred to 
as aeroisoclinicism, and its significance is that under these conditions 
the lift and center of pressure are substantially invariant with dynamic 
pressure. As pointed out in references 1 and 9^ there are certain dis- 
advantages attached to this condition; the lateral- control properties of 
such a 'Wing are net likely to be superior and the dynamic character- 
istics may well be inferior to those of a wing whicri is not operating 
at aeroisoclinic conditions. For the ordinary' wings represented in 
table 2, for instance, interpolation (on the reciprocals of the dynamic 
pressure at divergence) indicates that a wing with about 12° sweepback 
would be substantially aeroisoclinic at subsonic speeds; at supersonic 
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speeds the unswept wing is aeroisoclinic. That the unswept wing at 
supersonic speeds and probably also the wing with 12° sweepback at sub- 
sonic speeds are subject to large losses of lateral control can be seen 
from table 2. 

Some of the compounded wings represented in table 2 appear to be 
close to an over-all type of aeroisoclinicism because their dynamic 
pressure at divergence is very hi^. However, as a result of the more 
complicated geometry and the consequently more complicated structural 
properties of the compounded wings the fact that the dynamic pressure 
at divergence is approaching infinity is no longer a certain indication 
that the wing is approaching an aeroisoclinic condition. This subject 
will be discussed in some detail in the next section. 

In contrast to the compounded wings with very hi^ values of qjj, 

some of the compounded wings have dynamic pressures at divergence which 
are sufficiently low to be of concern, particularly plan forms of 
wings h, 5^ and 9* These plan forms, and probably plan form 3 also, 
must therefore be considered to be impractical from the aeroelastic 
point of view. 

In contrast to the dynamic pressure at divergence the dynamic pres- 
sure at aileron reversal is almost always of physical significancej 
unswept and particularly sweptback wings designed for high-speed flight 
are usually designed with resistance to reversal as one of the major 
structural design requirements. The dynamic pressure at reversal also 
serves as an index for the aeroelastic effects on the lateral control 
of a wing, but in itself it is only a crude index; for instance, althou^ 
the wing represented in figure 17 has a higher subsonic reversal speed 
than the one represented in figure l8, it has much less control power in 
the dynamic pressure range of primary concern (below q = 6) . As men- 
tioned previously this phenomenon may be predicted qualitatively from 
the ratio of the dynamic pressure at reversal to that at divergence. 

If, however, more complete information concerning the dependence of the 
rolling power and the maneuverability on the dynamic pressure is avail- 
able, the dynamic pressixre at reversal loses most of its significance. 

In view of the foregoing considerations no quantitative deductions 
should be made from the values of qj^ given in table 2. On the other 

hand, one conclusion may be drawn from them: Whereas the value of q^j 

given in table 2 vary from -13^ to +«, the values of q_ vary from 11.7 

£\ 

to 20.3 in the subsonic case and from 5-11 to 11.07 in the supersonic 
case; these numbers indicate that although the aeroelastic effects on 
the aerodynamic properties associated with level flight can be changed 
radically by a suitable compounding of the plan form the aerodynamic 
properties associated with rolling can be varied only within certain 
limits. That this is true for the process of balancing the effects of 
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bending and torsion deformations in general has been noted in refer- 
ences 1, 2, and 9- However, this conclusion must not be taken too 
literally; the limits within which the lateral control power can be 
varied by compounding are sufficiently wide to permit the selection of 
a satisfactory configuration in many cases. 


Extension of the Calculated Results to Other Plan Forms 

The process of compounding plan forms gives rise to two new geo- 
metrical parameters - the angle of sweep of the outer part of the wing 
and the position of the wing break - in addition to the three param- 
eters which define the geometry of the more conventional plan forms, 
namely, the angle of sweep, the taper ratio, and the aspect ratio. 
Although the 9 compo\mded plan forms considered in this paper are typi- 
cal of such plan forms, they fall short of the minimum number required 
to represent adequately all such plan forms that may be of interest. 

An attempt is made in this section to deduce the static aeroelastic 
characteristics of some related plan forms. 

The values of the parameters Cj /C t , Aa, C7 , and (— ) 

at a value of q = 3 are plotted in figure 23 as functions of the span- 
wise position of the break and as functions of the angle of sweepback of 
the outer panel in figure 2h. Also plotted are the lowest and second 
lowest (in absolute magnitude) values of q^. 

The series represented in figure 23(a) consists of plan forms which 
vary from a sweptback wing (y*^ = o) through a range of M-wings with 

Aj^ = -Aq = and varying positions of the break to a sweptforward 

wing (y*B = l-O) . At q = 3*0 the shift in aerodynamic center Ai. is 

seen to be 0 at both subsonic and supersonic speeds for the wing with 
y*B = 0 . 55 ^ and the loss in the lateral control powder at subsonic and 

supersonic speeds is zero for the wings with y*g =0.68 and 0.84, 
respectively. For all these wings Cj^ is positive, but for the wings 
with y*g = 0.55 and y*g =0.68 it is sufficiently high to be of no 
concern. ^The second- lowest value of q^ for subsonic speeds and for 
small values of y*^ is too large to be represented in the figure; it 
decreases from -80 to -120 as y*g increases from 0 to 0.3.) The rate 

of roll or lateral maneuverability is affected only slightly by a change 
in y*-gj as a result of the fact that changes in the rolling power are 

accompanied by almost equal changes in the damping roll. 
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The series represented in fignre 23(13) consists of plan forms which 
vary from a sweptforward wing (y*^ = o) through a range of W-wings with 

= -Aq = 45 ° and varying positions of the break to a swepthack wing 
(y *3 = 1 . 0 ). When y*g is about 0.6 the shift in aerodynamic center 

and the loss or gain in the lateral control power are almost zero at 
q = 3-0. However, the speed required to diverge this wing is relatively 
low. For wings with y*g less than 0.6 there is a gain in the lateral- 

control power due to aeroelastic action, but is even lower; for 

wings with y*g greater than 0.6 the dynamic pressure required to 

diverge the wings is higher, but there is some loss in control power. 
Again, the rate of steady roll is affected only slightly by a variation 
of y*g. 

The series represented in figure 23(c) consists of plan forms which 
vary from an unswept wing (y*^ = O) through a range of A-wings with 

A^ = 4-5°, Aq = 0 and varying positions of the break to an ordinary 
sweptback wing (y*^ = l-O)-. The aerodynamic -center shift is very small 
when y*g is less than 0.5, but this result is due to the fact that 

one-half or more of the wings is imswept. There is some loss in lateral- 
control power for all of these wings but the loss is very small when y*g 
is less than 0.3- The dynamic pressure at divergence is positive for 
most of the wings represented in figure 23(c) but sufficiently large to 
be of no concern. The rate of roll is substantially constant for y*g 

less than 0.4 and does not vary much for greater values of y*B* 

A series of A- type wings with A^^ = 45°, y*g =0.3 and varying Aq 
is represented in figure 24(a); when Aq is 0, they reduce to a A-wing 
with Ai = 45 °; when Aq is positive they are Intenaediate between a 
A-wing and an ordinary sweptback wing; and when Aq is negative, they 
are intermediate between a A-wing and a W-wing. At subsonic speeds the 
wing with Aq = -10° has no shift in aerodynamic center nor loss in 
lateral control; at supersonic speeds the wing with Aq = 0 has no 
shift in aerodynamic center, and the wing with Aq = -25° has no loss 

in control. The divergence speed of all three wings is probably suffi- 
ciently hi^ not to be of concern. The rolling speed does not vary much 
between the various wings represented in this figure. 

The wings represented in figure 24(b) differ from those represented 
in figirre 24(a) only in that their break is at 70 percent rather than 
30 percent of the semispan. The condition of zero aerodynamic -center 
shift can be achieved only at subsonic speeds in this series (with 
Aq = -22°) and the condition of no loss of lateral control power, not 
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at all. The rate of roll is about the sasie for all the wings repre- 
sented in figure 24(b) and only little lower than that of the wings 
represented in figure 23 and in figure 24(a). The divergence speed is 
sufficiently high for all the wings represented in figure 24(b) not to 
be of concern. 

A series of inverted A- type wings is represented in figure 24(c) j 
the wings differ from those represented in figure 24(b) only in that 
their inner parts are sweptforward rather than sweptback. The condi- 
tion of zero aerodynamic-center shift is not attained by any of these 
wings, and although their lateral- control power is more than adequate, 
their divergence speed is so low as to rule out most of these wings for 
practical purposes. 

The relation between the behavior of and the achievement of 

the aeroisoclinic condition can now be considered on the basis of fig- 
ures 23 and 24. In the case of ordinary swept wings the higher values 
of (corresponding to the higher modes) are much larger in absolute 

value than the lowest; a case where the two lowest values coincide in 
absolute magnitude does not appear to arise for most ordinary swept 
wings. Consequently, if for these ordinary swept wings ^ is plotted 

as a function of the angle of sweepback or, more generally, as a func- 
tion of the parameter k defined in references 1, 3^ 4, and 9 (or the 
parameter d/a of ref. l4) , which contains the stiffness ratio GJ/EJ, 
the aspect ratio and the moment arm e-j_ in addition to the angle of 

svfeepback, there is only one value of the parameter for which q-^ goes 
to infinity. 


For the compounded wings, however, the two lowest values of qjj 

frequently have the same absolute magnitude. For instance, in fig- 
ure 23 (a) the lowest value of is negative at subsonic speeds for 

all values of less than 0.59 and the second lowest value is posi- 

tive for values of y*g greater than about 0.35 and less than 0.59- At 
y *3 = 0.59 the two values of qjj coincide in absolute magnitude, and 
at values of y*g greater than 0.59 the lowest value of is posi- 

tive and the second lowest value is negative. Consequently, the lowest 
value of never approaches infinity; it merely jumps from negative 

to positve at y*-g = 0.59. The same phenomenon occurs at supersonic 
speeds at y*g = 0.71- The second lowest value of Ojj goes to infinity 
twice, at y*g = 0.35 and 0.75 at subsonic speeds and y*^ = 0*50 
and 0.87 at supersonic speeds. 
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For the family of wings represented in figure 23(a) the aerodynamic- 
center shift was noted to be 0 at y*g = 0.55 for both subsonic and 

supersonic speeds. The change in lift-curve slope due to aeroelastic 
action is 0 for y*g = 0.53 for subsonic speeds and y*g = 0.63 at 

supersonic speeds. The conditions of either zero aerodynamic- center 
shift or zero lift increase may be considered to define the over-all 
aeroisoclinic condition. Inspection of figures 23 and 2k indicates 
that the jump of the lowest value of from negative to positive 

tends to occur when the wing is close to an over-all aeroisoclinic con- 
dition, but more definite conclusions cannot be drawn. 

The values of Qjj shown in figure 24(a) vary with the angle of 

sweep of the outer part of the wing in much the same manner as they do 
for ordinary swept wings with the angle of sweep of the entire wing. 

The lowest value of goes to Infinity at = 12° and -2° at sub- 

sonic and supersonic speeds, respectively. These values of Aq are also 

the angles at which the change in the lift and the shift of the aero- 
dynamic center are zero, although at subsonic speeds the aerodynamic- 
center shift is also zero at Aq = -l8°. On the other hand, the values 

of qp shown in figures 24(b) and 24(c) vary with Aq in an entirely 
different manner. The lowest values of q^j never change sign; the 
second lowest values of Op go to infinity at Aq = 21° and Aq = -6° 

at subsonic and supersonic speeds, respectively, in figure 24(b) and 
at Aq = 33° and 0, respectively, in figure 24(c). However, these 
values of Aq have no significance insofar as the lift, aerodynamic- 
center shift, and lateral control power are concerned, as may be seen 
from figures 24(b) and 24(c). Consequently, any deductions concerning 
the aeroelastic phenomena of interest can be drawn only from the behav- 
ior of the lowest value (in absolute magnitude) of q^-j. 


The Optlniim Compounded Plan Form 

On the basis of the preceding discussion the problem of the selec- 
tion of an optimum compounded plan form can now be broached, the term 
optimum being used in the sense of most favorable aeroelastic character- 
istics at the least sacrifice in aerodynamic and structural performance. 
In view of the relatively small number of plan forms considered in this 
paper and inasmuch as no dynamic aeroelastic calculations have been made, 
the following discussion can shed light on only a few aspects of the 
problem. 

The desired static aeroelastic characteristics .- The desired static 
aeroelastic characteristics are, approximately in the order of their 
importance : 
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(1) The shift of the aerodynamic center should he small if it is 

forward, as it is in almost all cases of practical interest. 

(2) There should, he no loss in the rate of roll nor in the lateral 

control power. 

(3) There should he no appreciable change in the lift-curve slope. 

(4) The dynamic pressure required for divergence should he either 

negative or, if positive, at least 25 percent higher than 
the highest expected dynamic pressure in the given speed 
range . 

In the case of airplanes designed to fly occasionally at high supersonic 
Mach numbers ( say 2 or greater) these conditions should he satisfied as 
much as possible at those Mach numbers as well as at subsonic speeds. 

The selection of a plan form possessing some of these character- 
istics can now he discussed on the basis of a union effected between the 
three series of plan forms represented in figure 23 with those repre- 
sented in figure 24, a process which gives rise to plan forms of which 
the inner part is either swept bac2< or swept forward at an angle of 4-5°, 
but of which both the outer part and the position of the wing break are 
arbitrary. In considering all these combined results the fact must be 
kept in mind that they still apply to a quite restricted class of wings. 
The aspect ratio of all w'ings is 6, the taper ratio 0.5, and the angle 
of sweep 45 ° or -45°, the span of the outboard aileron is 50 percent of 
the wing span (unless another aileron configuration is specified) , and 
the structures are of a certain kind, namely the stiffness distributions 
vary substantially as dictated by constant- stress considerations (as 
outlined in ref. l) , and GJ/EI = 0.794. The conclusions reached in 
this section may not be valid for any wing -with taper ratio, aileron 
configuration ( in the case of conclusions concerning lateral- control 
properties) , stiffness ratio, and variation of the stiffnesses GJ 
and El (particularly near the wing tip) which differ greatly from the 
values used in the calculations described in this paper. 

The information pertinent to this discussion is summarized in fig- 
ures 25 and 26 . The plan forms represented in figures 25(a) and 26(a) 
have an inner part sweptback 45°, and those represented in figures 25 (b) 
and 26 (b) have an inner part which is sweptforward 45°. The location of 
the break and the angle of sweepback of the outer part are arbitrary and 
constitute the coordinates of these figures. The curves shows;! are the 
locr: of the points representing plan forms which have a zero aerodynamic- 
center shift at q = 3-0 (in the case of fig. 25 ) and no loss in lateral 
control at q = 3*0 (in the case of fig. 26). The aerodynamic-center 
shifts for all wings represented in figure 25 are negative, except that 
those represented by points on or between the two lines for zero shift 
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have a zero or a very small positive shift, respectively. The rolling- 
moment ratio is less than 1 in the region above and greater in the 
regions below the lines of no loss of lateral control shown in figure 26. 

A comparison of figure 25 with figure 26 indicates that a wing with 
the inner part swept forward ^ 4 - 5 ° must have the break at 72 percent span 
and the outer part swept back 6o° in order to have no shift in aero- 
dynamic center nor loss in lateral control at subsonic speeds. No wing 
with the inner part swept forward ^5° can satisfy both of those condi- 
tions at supersonic speeds. When the inner part of the wing is swept 
back, however, almost all the wings represented by the line for subsonic 
speeds in figure 26(a) should be satisfactory at subsonic speeds because 
they have no loss in lateral control and only a small forward shift of 
the aerodynamic center. Of this group of wings those on the lower part 
of the curve (y*g greater than 0-55 or Aq less than -4o°) should be 

satisfactory at supersonic speeds as well as subsonic speeds because 
they should have only a small rearward shift of the aerodynamic center 
and only a small loss in the lateral control. Probably the optimum 
wings among those considered here are the ones with y*-g = 0.58 to 0.6o 

and corresponding values of Aq from - 40 ° to - 45 °. 

The mmiber and range of geometric and structural variables covered 
by the calculations described in this paper is insufficient to permit 
any generalization of the figures cited in the preceding paragraph. In 
a given case sufficient calculations should be made to permit the prepa- 
ration of charts similar to those of figures 25 and 26 for several values 
of the sweep angle of the inner part of the wing and, unless they are 
decided ’upon beforehand, for several aspect and taper ratios. Also, if 
the stiffness-ratio GJ/eI can be varied without increasing the struc- 
tural weight, several values of the ratio should be considered. Further- 
more, inasmuch as the simultaneous achievement of zero aerodynamic-center 
shift and zero lateral- control loss at both subsonic and supersonic 
speeds in an aerodynamlcally acceptable wing is \inlikely, it would be 
desirable to plot on charts of the type represented by figiares 25 and 26 
lines of constant aerodynamic-center shift and lines of constant lateral- 
control loss, respectively, in addition to lines of zero shift and zero 
loss in order to facilitate the selection of an optimum compromise plan 
form. (The number of plan forms for which calculations have been made 
is insufficient to permit the plotting of such contour lines on figs. 25 
and 26.) This procedure implies a great number of calculations, despite 
the fact that many plan forms with obviously undesirable aerodynamic 
characteristics can be eliminated from consideration beforehand, as will 
be pointed out in a subsequent section. 

The no-shift and no-loss requirements can, of course, also be satis- 
fied simultaneously by choosing a wing with zero aerodynamic-center shift 
and equipping it with an all -movable -wing tip. This procedure has the 
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advantage of providing greatly increased rigid-wing rolling performance, 
particularly at supersonic speeds, in addition to the decreased loss in 
this performance as a result of aeroelastic action. However, this 
advantage is offset to a large extent hy the mechanical and flutter 
problems which beset such a configuration. Only if the wing under con- 
sideration has its break far outboard, say 85 percent or more of the 
semispan, and the lateral control provided by rotating the portion of 
the wing outboard of the break is sufficient is this configuration 
likely to be practical. 

The desired dynamic aeroelastic characteristics .- The desired 
dynamic aeroelastic characteristics are, substantially, 

( 1 ) That the inertia effects on the aerodynamic center and the 

lateral control should be small or in such a direction as 

to oppose any adverse static aeroelastic effects 

( 2 ) That the flutter speed be higher thsui any expected flying 

speed at all altitudes 

( 3 ) That the dynamic response of the wing to atmospheric excita- 

tion give rise to no excessive stresses 

The inertia effects on the aerodynamic -center shift will be small 
if the wing weight represents a small fraction of the aiiplane weight, 
and the inertia effects on the lateral-control power are not likely to 
be important because the rate of steady roll (which is independent of 
inertia effects) is usiially considered to be more important than the 
control power (which is an index of the attainable rolling acceleration) . 

The flutter and dynamic-response characteristics of wings designed 
on the basis of a balance of torsion and bending deformations may well 
be inferior to that of an ordinary swept wing because they are likely 
to be more flexible. Careful flutter studies must therefore be made in 
each case. When necessary, however, it may be possible in some cases 
to raise the flutter speed at relatively low weight penalty by taking 
advantage of the large moment arms available for mass balancing. A 
high wing- flutter speed (relative to the highest flying speed at the 
given altitude) usually implies satisfactory dynamic -response charac- 
teristics, provided the airplane as a whole is stable. 

The desired structural characteristics .- For the purposes of this 
discussion the desirable structural characteristics are that the weight 
of the structure of a compounded wing be no higher than that of a com- 
parable ordinary swept wing. The break requires locally a certain 
amount of material not needed in an ordinary wing, and the large torques 
near the root of an M- or W-wing require additional torsion-resisting 
material. However, the bending moments are much smaller than for an 
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ordinary wing, and the saving in flexure-resistant material may he so 
large as to compensate for the aforementioned increases in structural 
material. Thus the weight of a compounded wing may he little, if any, 
higher than that of an ordinary wing. 

The desired aerodynamic characteristics .- Within the scope of this 
discussion the desired aerodynamic characteristics are that the drag 
and the stability characteristics of the compoiuided wing he no worse 
and better, respectively, than those of a comparable swepthack wing. 

A few tests have shown that the stability characteristics of M- and 
W-wings can, indeed, be superior to those of ordinary swept wings. The 
drag of a compounded wing is likely to be higher than that of an ordi- 
nary swept wing, as has been shown by tests at subsonic, transonic, and 
supersonic speeds (see ref. 15, for instance) and as may be inferred 
from calculations for ordinary swept wings at supersonic speeds (see 
refs. l6 and IT, for instance). The whole question of whether to use 
M- and W-wings thus hinges primarily on the problem of whether the 
saving in structural wei^t afforded by these configurations in achieving 
the desired stability and control characteristics is worth the drag 
penalty. 

The additional drag of a compounded wing as compared to an ordinary 
sweptback wing arises from three sources: the fact that part of the 

wing may have a relatively low sweep angle, the fact that in the case 
of an M-wing the inner part of the wing may be sweptforward (giving rise 
to fuselage interference drag) , and the fact that in the case of a W-wing 
the break itself is the source of a certain amount of interference drag. 

When the results of static aeroelastic analyses of a variety of 
plan forms are summarized on charts similar to figures 25 and 26 some 
of the plan forms brought to light by these charts can be eliminated 
because they are likely to be subject to one or more of the afore- 
mentioned types of drag. In figures 25 and 26, for example, the plan 
forms represented by points within the wedge-shaped regions labeled 
"aerodynamically undesirable" have too much of their area swept at too 
low an angle to compete with a completely swept wing. (These regions 
are based on qualitative estimates and are shown primarily for illus- 
trative puip)oses.) 

Similar reasoning may be employed in connection with the inter- 
ference drag. In figures 25(a) and 26(a), the plan forms represented 
by points above the wedge-shaped regions are likely to have higher drag 
than do those represented by points below this region, because the 
difference in sweep angle between the inner and outer portions of the 
wing is much greater for the former than for the latter. The reverse 
is true for figures 25(b) and 26(b), Also, the plan forms represented 
by points in figures 25(a) and 26 (a) are likely to have less drag than 
those represented by points in figures 25 (b) and 26 (b) because the 
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interference drag caused 'ey the break of a W-wing is likely to be less 
than that caused by the root of an M-wlng as a result of the smaller 
chords involved, particularly if the break of the W-wing is near the 
xip. 


Finally, experimental or theoretical drag studies must be used to 
decide which, if any, of the remaining plan forms may be suitable for 
any given purpose, the theoretical studies being useful primarily for 
supersonic speeds. If there are any satisfactory compounded plan forms, 
these studies should be followed by further studies aimed at reducing 
the drag of these wings. For instance, there is a possibility that the 
interference drag at subsonic speeds may be reduced by using fences. 
Also, the drag caused by the region of the wing which has a relatively 
low sweep angle can be reduced by resorting to thinner airfoil sections 
in that region; the resulting w’eight penalty should be very small, 
because that region is likely to contain only a small part of the wing 
area and be in a region w'here the stresses are relatively low. Thus, 
when large adverse static aeroelastic effects are anticipated, as for 
wings with low wing loading designed for low load factors and intended 
for high-speed low-altitude flight, compounded plan forms may well 
constitute the best solution. 

CONCLUDING REMARKS 

Calculations have been made of the static aeroelastic character- 
istics of nine M-, ¥-, and A-wings by using the best available aero- 
dynamic and structural information. Although the number of plan forms 
considered is too small and the calculations too specific to permit of 
quantitative conclusions which are generally applicable, certain quali- 
tative conclusions have been draw^n. The question of the plan form with 
the optimum static aeroelastic characteristics has teen discussed on 
the basis of these conclusions. In general, there is reason to believe 
that by suitable compounding plan forms can be obtained which are 
superior aeroelastically and structurally to ordinary swept wings. 


Langley Aeronautical Laboratory, 

Rational Advisory Committee for Aeronauxics, 
Langley Field, Va. 
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APEENBIX 

METHOD OF CALCULATING STATIC AEROELASTIC PHENOMENA 
OF M-, W-, AND A-WINGS 
Symaietrical Flight Conditions 


The aerodynamic Influence coefficients .- In keeping with the assump- 
tions concerning the aerodynamic properties of the wing, the lift on sec- 
tions parallel to the air stream is given hy 



( 1 ) 


where is an aerodynamic influence-coefficient matrix. A method of 

calculating such a matrix from known rigid-wing additional lift distri- 
butions is given in reference 5- The matrix obtained in this manner is 





+ 



( 2 ) 


where 



is a unit matrix. 


is a 


square matrix defined by 


1 0 ... 0 

10 . . .0 

[ir] = (3) 


1 0 . . . 0 
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and 



a column matrix, each element of which is unity. 


The row 


matrix 



serves to integrate the lift distribution and is based 


on Simpson's rule with a modification which assumes that the lift goes 
to zero with Infinite slope at the wing tip. (See ref. 3 *) The param- 
eter is given in reference 5 in terms of the plan-foim parameter 


F = 


2it 


a 


cos A 


by the relation 






F./l + 




+ 6 


ih) 


For M-, and A-wings this method requires some modification. 

For the purpose at hand the expression for [qJ (eq. (2)) can be 
written as 



In reference 5 tbe factor K (which may be considered to be the ratio 
of the lift coefficient due to a unit symmetrical twist to the dimen- 
sionless lateral center of pressure of the additional lift distribution) 
and the related factor K' used for anti symmetrical lift conditions are 
shown to be 1 according to lifting-line theory, and the method of calcu- 
lating approximate aerodynamic influence coefficients presented in refer- 
ence 5 is based on the assumption that these factors are 1. The calcu- 
lations made in reference 6 have shown, however, that these factors are 
not 1 for M-, W-, and A-wings. The values of K and K' are therefore 
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obtained from the spanwise lift distributions and associated aero- 
dynamic parameters presented in reference 6. They are introduced in 

the matrix [q^ (respectively the matrix w in the antisymmetric 

case) in such a way (see eq. (2)) as to yield the correct spanwise lift 
distribution for any angle-of-attack distribution which consists of a 
linear superposition of a constant angle of attack and a linear twist 
in the symmetric case, and of a linear twist and a 50 -percent- semi span 
outboard-aileron deflection in the antisymmetric case. Specifically, 



Also, inasmuch as the plan-form parameter F is a function of the 
cosine of the sweep angle, there is some question as to which value 
should be chosen in the case of a A-wing or of an M- or W-wlng with 
angles of sweep in the inner and outer parts of the wing which differ 
in absolute magnitude. For the calculations made by the method of the 
present paper an average value of cos A deduced by 

c 

cos A = cos A — dy* 

^0 c 


has been used. 

This procedure is believed to furnish results with sufficient 
accuracy for the purpose intended, because the values of the lift dis- 
tributions calculated by the matrices and are not very 

sensitive to the value of F. Since, when K and K' are 1, the 
matrices [q^ and reduce to those presented in reference p, 

which are valid for all angle-of -attack distributions, there is reason 

to believe that the matrices and used in the present paper 

yield lift distributions for angle-of-attack distributions other than 
the aforementioned ones with sufficient accuracy for the purpose 
Intended. A few calculations by means of the method used in refer- 
ence 6 and of the matrices of the present paper for parabolic symmetric 
and antisymmetric twists have yielded results in excellent agreement 
with each other. 
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If strip theory is used, as has been done in the calculations for 
supersonic flow, the aerodynamic influence-coefficient matrix is 


H - Lt] 

The aeroelastic equation .- The running load on each part of the 
wing can be written as 

{z} = (7) 

and the running torque in planes parallel to the air stream for 
uncambered wing sections is 



( 8 ) 


The subscripts i and o used on the matrices of equations (7) and (8) 
in ohe following analysis refer, respectively, to the inner and outer 
parts of the wing. 


The single and double integrations required to obtain the accumu- 
lated torques and bending moments from the running torques and loads 

are performed by means of integrating matrices H and , respec- 

tively. These matrices are based on Simpson's rule and are similar to 
nhose described in reference 3* When a modification is made at the 
wing tip which takes into account the infinite slope of the spanwise 
lift distribution at the wing tip at subsonic speeds, this fact is 
described by adding a prime mark to the matrices. For the sake of 

definiteness the matrices [l [] and [ll G appropriate for subsonic 
speeds will be used in the following derivation. 


The integrating matrices used in the calculations described in the 
present paper are given in table 3* They pertain to the stations used 
in the calculations, namely: 


For ordinary wings: y* = 0.1, 

For wings with y*g =0.3: y* 

0.86 


0.25, O.it, 0.55, 0.7, 0.85 
= 0.1, 0.2, 0.3, 0.4i^, 0.58, 0.72, 


For wings with y*^ 
0.75, 0.875 
For wings with y*^ 


= 0 . 5 : 

y* = 0 . 1 , 

0.2333 

= 0 . 7 : 

y* = 0 . 1 , 

0 . 25 , ( 


0 . 3667 , 0 . 5 , 0 . 625 , 
0 . 55 , 0 . 7 , 0 . 8 , 0.9 
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In the calculations of the aerodynamic influence coefficients, stations 
were taken at y* = 0, 0.05j and 0.1 on the part of the wing covered 

by the fuselage for all wings. The matrices [_lfj 
shown in table 3j they are the same as the matrices |_-;J . 

for the wings with y*g = 0.3j except that when used as matrices LifJ 
and they refer to the stations y* = 0, 0.05j and 0.1 rather 

than the values given in table 3j namely y* = 0.1, 0.2, and 0-3- 


The accumulated bending moment M about axes parallel to the air 
stream can then be obtained directly by using the matrix CnO to 

perform the indicated double Integration. Similarly, at a given sec- 
tion parallel to the air stream the accumulated torque about an axis 
which passes through the shear center of that section and is perpen- 
dicular to the plane of symmetry can be found by performing a single 

integration by means of the [!■] matrix of the running torque and 

then subtracting from this result the product of the accumulated bending 
moment at the section and the tangent of the angle of sweepback outboard 
of the section. 


Thus 



( 9 ) 


( 10 ) 


- tan . (ll) 

- tan (12) 

The prime marks on and . indicate that the moments carried 

across the wing break from the outer part of the wing are not included. 

In these equations the station at the wing break may be considered 
to be the last one on the inner part of the wing or the first one on the 
outer part of the wing or both. In the following derivation the last- 
named alternative is implied, except where specified otherwise. 
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In order to find t;he oending and twisting deforn^tions along the 
wing, the accumulated bending ir.omen~ about axes perpendicular to 

the elastic axis and the accijinulated torque about the elastic axis 

must be known. These new moments at any station may be obtained by 
means of the transformation 


cos A -sin A I M 


sin A cos A It 


which yields 




{-a}, = ^ I 


{“a} = - Ao I 

c 

^’a} ■ Al I <=°s 

{’'a}o * '"o I AoCl]^^ 


When the moments carried across the wing break from the outer part of 
the wing are added to the moments [“'a], and ^ _ the following 

expressions are obtained: 

{Ma}. = {m'a}. ^°s(^i - ^'^o)||1b|!J?’^a}^ ' si^(Ai - Ao)||1b||.{ta}^ 

0} (18) 

*1 1 O O 
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sin(A^ - Aq) + cos(Ai 



(19) 


where -the raatrix 



is a rectangular matrix defined as 



1 0 . . . 0 

10 . . . 0 


1 O . . . 0 


( 20 ) 


in which the numher of rows is equal to the nimiber of stations on the 
inner part of the wing and the numher of columns is equal to the numher 
of stations on the outer part of the wing. The diagonal matrix 


is composed of the dimensionless moment aimis of the normal shear at the 
wing break about the stations on the inner part of the wing. The temn 
of equation (l8) in which it occurs represents the contribution of the 
vertical shear at the wing break to the bending moment at sections on the 
inner part of the wing. 


The preceding equations may be combined as follows: 





( 21 ) 


( 22 ) 
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the matrices and jjr^ (with subscripts) being defined as 


nj."] = [if] - d-; ^ sin [i] 

LJii cos A-l i 13/2 


S2_c 


( 25 ) 


Mlio - M.BG, 


^ ^ =ln Aillsli [l JJ^] 


C 3 , 

b/2 


Hi] = ^2 [ii'] - 2l- sinAo[l’] 

LJoo cos Aq ^ -'o ^ b/2 °L?r 1 


(26) 


(27) 


prl = d^ cos Aiffl 

‘--'il -b /2 -^-^i 


e-. c 


(28) 




(29) 


Heo = ^o^-A„ClO„[gi| 


( 30 ) 


According to simple beam theory, the equations of equilibria of 
the deformed wing are 


GJ 



= Ta 


( 31 ) 
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for twisting about the elastic axis and 


El 



(32) 


for bending along the elastic axis. These equations may be integrated 
to obtain qp and T with the matrix [X| ' ' given in table 3- This 
integrating matrix is (as explained in ref. 3) the double transpose of 
[l] if the stations taken along the wing span are equally spaced. 

The wing root is considered to be the wing section perpendicular to 
the elastic sLxis which passes through the intersection of the elastic 
axis and the side of the fuselage. Unlike the case of an unswept wing 
there are rotations of this section when the wing is subjected to bending 
moments and torques. These rotations have the nature of rigid-body rota- 
tions imparted to the rest of the wing and are caused by the flexibility 
of the root triangle and of the can^- through bay inside the fuselage. 
They can be calculated by analyzing these components in detail and can 
be expressed in terms of the four dimensionless flexibility constants 
defined in reference 3 j 


\ 




We/(GJ)r 

w 

'^e /( ) r 


Wg/(EI)^ 


We /(SI) I 


(33) 

(35) 

(36) 


■Wg being defined as in figure 1. They must then be added to the values 
of qp and r obtained from equations (3l) and (32). As pointed out in 
reference 3 } the values of (GJ)^. and (El)j. serve only for reference 

purposes j so that their exact values are immaterial. The values obtained 
by extrapolating curves of GJ and El plotted along the span to the 
root station are probably the most convenient ones to use. 
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The resulting expressions for the wing deformations are 




2 


(GJ);^ cos Aj^ 


m Al _ Pi 

®i Lgj li 




Wg cos 

I h 

dl 5 


%KlWj 


(37) 


d u 

n T" 


fcp> = :i-i — [i] ' 

(GJ)^ CCS Aq ‘--^o 


(GJ). 


GJ 




(38) 




-I 


(EI)j, cos A^ 


‘U' 


Hi’ 


(SI). 


El 


Wg cos A 




^ 2 


Wg cos A^ 






£«a}. 


(39) 


and 


© - 


.<3-0 2 


H, 


(El)j. cos Ag^-^o L 


(GJ), 


GJ 


{“a}. 




where the values of cp and F on the outer panel are measured rela- 
tive to the station at the wing break rather than to the station at the 

wing root. The matrix Is a square matrix w'hlch has as many 
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columns as these stations on the inner part of the wing and is defined 

ty 



0 0 0 . . . 0 

10 0 . . - 0 

10 0 . . . 0 


10 0 . . . 0 


Rigid-body- rotation constants could also be introduced to take 
into account the effect of local distortions in the vicinity of the 
wing break on the deformations of the outer part of the wing. No such 
constants have, however, been used in the calculations described in 
this paper because no simple method of calculating them was available , 
and if they had been calculated for a specific case, the results would 
not necessarily have applied to other cases. Also, inasmuch as they 
affect only the deformations of the outer part of the wing, whereas 
those considered in the preceding paragraphs affect the entire wing 
(and even those do not have a large effect on the wing deformations 
except for wings of low aspect ratio) , there is good reason to believe ‘ 
that the rigid-body-deformation constants appropriate to the wing break 
can be neglected. 

iihe angle of attack due to structural deformations ag is related 
to q) and P (for small angles) by 

. = cos - sin (^l) 


= |11b||o{;^s}. + cosAo{?}^ - sinAo{r}^ 


where the rectangular matrix |lg|| defined by 


0 

0 


0 1 
0 1 


I^bII = 


0 1 


(1^2) 


{h3) 
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has as many rows as there are stations on the outer part of the wing 
and as many columns as there are stations on the inner part of the wing. 

The substitution of equations (37) and (39) into equation (4l), 
and of equations (38) and (itO) into equation (42) yields 




(44) 






o o 


(45) 


where the matrices H are defined hy 


[Kt]^ = 


(GJ), 


GJ 




w_ cos / (GJ)_ 

+ — : f 0^ - T tan AiQp, 


h/2 (El). 




(46) 


n" “1 (GJ) I — I 1 

Hi ■ (ST 


(El). 


El 


Wp cos A,- / (GJ)t;. 

® ^ ^ ^ tan A. 


b/£ 


-K- 


(EI), 


i%)&0 


(47) 




(48) 


r -1 

^ (5)7 


tan A^Cl] 


(El) 
El 


(49) 
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The comhination of equations (2l), (22), (23), (24), (44), and (45) 
then gives 


y 

_ .=4)' 

HiiWIio 

[al 


(C-J)r 


L>4oj 


(50) 


where the matrices [b^ and j|B|[ (with subscripts) are defined as 

[B]ii - (51) 


|B|l„i - ilB||„CB],^ 

Hoo - INlJl^lio * HoHoo - WoH 


oo 


(52) 

(53) 

(54) 


Equation (50) may be written as 



(55) 


The matrix jjB^ is defined as 




( 56 ) 
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In confining the four suhmatrlces the order of the resulting matrix 
may be reduced by one because the station at the break is represented 
twice - by two rows and two columns. The combination is effected by 

omitting the first rows of the submatrices IsIq-; &] oo and by 

adding the first column of the submatrices H^llj[^o oo 

last columns of submatrices and respectively. 

Substitution of equation (l) into equation (55) yields the desired 
aeroelastic equation for symmetric flight 



(57) 


where the dimensionless parameter q is defined by 



and the aeroelastic matrix [a] is defined by 

H =[b][q] 


(58) 


(59) 


Solution of aeroelastic equation, for dynamic pressure at diver- 
gence . - The condition for divergence is that the elements of Be 

finite when the geometric angle-of-attack -S zero along the 

entire span; therefore, the value of the parameter q at divergence 
is the lowest real positive value of q^j which satisfies the equation 



(60) 


If the lowest root is real and distinct, it can be computed by simple 
iteration. Often, however, the dominant roots of the matrix [a] for 
M- and W-wlngs are not very well distinct and the simple iteration 
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procedure does not converge. If the matrix [a^ does have two domi- 
nant roots distinct from the others, the iteration procedure converges 
rapidly to the following relation between successive iteration columns: 


where are the dominant characteristic values and 

{<^} is the trial column for -^a^. This equation represents as many 
linear algebraic equations as there are rows or columns in the matrix [aJ . 
Each equation involves the two unknowns ^ 




and any two of the equations can be solved for these unknowns, and, hence, 
for and • (See also ref. l8-) 

If the' simple Iteration of equation (6o) yields a small negative 
value of there is the possibility that for these compounded con- 

figurations the next larger values (in the absolute sense) may be posi- 
tive and still low enough to be of concera. The next higher value of 

qjj may be found in the following manner. The modal column 
obtained by the simple iteration of eqixation (6o) is orthogonal to a 

(2) 

modal row [_pj ' ' which corresponds to the second mode of divergence 
and satisfies the relation 




The orthogonal relationship between •C“"}w 
expressed by 


( 1 ) 


and 




(1) = 0 


or 




( 
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where the matrix 




is defined by 


0 0 
0 0 


>] = Lil - 



0 



( 1 ) 


0 0 




(i)L ( 1 ) 


’n-l 




=n 


0 0 ... 0 


1 


The substitution of 



for 


yields 


LfJ 


( 2 ) 





The first node of divergence has been eliminated in this equation; hence 
~ (2^ 

tne value of may oe obtained by simple iteration. Tne correct 

value of i^s-y also be obtained by the iteration of HM with 

a column matrix, although the modal column obtained is spurious in the 
sense that it will not satisfy equation (6o). (See also ref. l8.) 

Solution of the aercelastic equation at subcritlcal conditions . - 
E quation (57) may be rearranged to read 


because 



(61) 


(62) 


The total angle of attaclc can then be found for any type of geometrical 
angle of attack by solving the simultaneous eqixatlons represented by 
matrix equation (6l), using the given geometrical angle-of-attack distri 
buttons as the "knowns" on the right side of equation (6l). Then the 
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corresponding lift distribution is obtained by premultiplying the angle- 
of -attack distribution by the matrix H- For instance, the additional 
lift distribution for the flexible wing is obtained by setting equal 

to unity along the entire semispan in equation (6l) and then premulti- 
plying the resulting column by namely 



( 63 ) 


Within the limitations of the assumptions in the derivation of the 
matrix (see ref. 5) j "the lift distribution on that part of the 

wing covered by the fuselage is proportional to the rigid-wing lift 
distribution in that region and is determined in magnitude by the angle 
of attack at the wing rootj specifically. 



(64) 


Integration of the lift distribution represented by equations ( 63 ) 
and (64) yields the lift coefficient. 





cc^ 


c^C; 


L 


S) 




(65) 
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and the bending-moment coefficient 


Cg - 


qSb 




CC: 


f c„C 


r'-L, 


°0 


+ d. 






cc • 




'«0 


> + 
i 


* (% * DtllJoJ 


r 


cc ■ 


'^rCr^ 


o> 

(66) 


where the row matrices [jiJ first rows of the 

integrating matrices jjl^ and respectively. 

SiiT-ilarly the pitching-moment coefficient is 


C’Tl — 


_ 


qSc 




CCl 


^r^L, 




+ d. 




X 


cc : 


CrCr 


'«0 


> + 
i 




cc ] 


^r^L. 




( 67 ) 


the parameter x being the streamwise distance of the local aerodynamic 
center rearward of an unswept reference line through the quarter-chord 
point of the mean aerodynamic chord. 

The lateral center of pressure ^ and the position of the wing 
aerodynamic center a may be calculated from the lift and moment 
coefficients given by equating ( 65 ), (66), and ( 67 ) as fellows: 
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(68) 


a = — - — 

4 Cl 


(69) 


Anti symmetrical Flight Condition - Damping in Roll 

The aerodynamic influence coefficients .- Equation (l) holds for 
antisymmetric lift distributions provided the aerodynamic influence- 

coefficient matrix is replaced hy a matrix M- A method for 

calculating this matrix is given in reference 5 but, as is tznie for the 
symmetric lift distributions, certain modifications have»to be made in 
the case of M-, and A-wings. The modified matrix is 
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and where the constant K' obtained from the rigid-wing lift distri- 
butions due to aileron deflection and due to antisymmetrical linear 
twist is 



The aeroelastic equation .- Equations (57) and (6l) apply to anti- 
symnetrical loadings provided the aeroelastic matrix [A^ is replaced 

by a matrix M defined by 


[Aa] = [B][0a] 




Solution of •Qhe aeroelastic equation .- The unit antisymmetrical 
linear lift distribution for the flexible wing is obtained by setting 
ttg equal to y* in the antisymmetrical equivalent of equation (6l) 

and premultiplying the resulting column -{^c^ by > 3-s follows 




w 



(75) 


Equation (64), which gives the symmetrical lift on that part of the wing 
span covered by the fuselage, also applies to antisymmetrical lift and is 



cc ■ 




°0 
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the parameters l-r, T-t^, and now pertaining to the anti- 

0 c^Cl^ 

symmetrical case. 


The damping- in-roll coefficient can he obtained by integrating 
this spanwise lift distribution, so that 



( 76 ) 


Anti symmetrical Flight Condition - Wing Loading 
Due to Aileron Deflection 

The aeroelastic equation .- The manning load on a flexible wing with 
aileron deflected is 


'C^} = 




( 77 ) 


where the superscripts (°-s) and (S) refer to the lift distributions 
due to wing deformation and unit aileron deflection, respectively. The 
corresponding running torque is 



( 78 ) 


V 


U8 
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The subscripts 1 and o used in the following analysis on the matrices 
of equations (T6) and (77) refer, again, to the inner and outer part of 
the wing, respectively. 

The use of the two preceding expressions in equations (l^i) to (2^0 * 

yields 


r' 
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where the matrices jl and [rj (with subscripts) are the same as 


the corresponding matrixes M and U] , except that 


eic 


must be 


replaced everywhere by the diagonal matrix 

2 


eoc 


; thus , 


H [ll] + d. ^ sin A, [l] 

‘--'ii cosAi' 'i ^ b/2 ^^"^i 


egC 


(83) 


do 






(84) 


&] = — ^ ^ ^ ^o&O (85) 

^-^oo cos Aq ^ ^ b/2 o' ^o[cj. j 


Hii = ^ 


eoc 


c^ 


(86) 


2 sin(Ai - -do) 


■"“lo = "”cos A, 


ll.&iD, 


d, cos Ail[lBll.[l-]^g^ 


(87) 
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Substitution of equations (79) to (82) into equations (^4) and (45) 
yields 



Hoc - + &0<,Hoo - 


(93) 
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Equation ( 89 ) may then be rewritten as 


^ Vv - (GJ) 




- B 


cc ■ 


(S) 


w 


(94) 


where 


B = 


^iiPIlio 

b11 -0 

"01 *- -*oo 


(95) 


except that, again, in the processes of constructing the matrix [Bj 
from its submatrices, it is reduced in order in the same manner as the 
matrix [^bJ . 

The aeroelastic equation for aileron deflection is then obtained by 
substituting equations (74) and (75) in equation (94), so that 


I “ 0 ^ 


(5) 


( 96 ) 


Solution of the aeroelastic equation .- Equation ( 96 ) can be solved 
in the same manner as equation C^l). Once it has been solved for a given 
set of knowns defined by the right side of equation ( 96 ), the rolling- 
moment coefficient due to a unit aileron deflection can be obtained from 


2M^ 


(5) 


' ^5 qSb 


= C 



1 

2 








i 


+ 



(97) 
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where 




( 98 ) 


The rate of roll per iinit aileron deflection is then given by 



( 99 ) 


If C 2 g is for -unit effective aileron deflection ^unit y instead 
of unit actual aileron deflection then 

(i), , - 5 

o=x ‘•p 

The condition for aileron reversal is that the rolling-moaient coef- 
ficient C2g be zero. Combining the row integrating matrices of equa- 
tion (97) into one row matrix, and setting the resulting expression 
equal to zero yields 




= 0 


( 101 ) 



where , as in the case of the matrices [b] and , uhe order of the 

first row matrix on the right side of equation (l02) is of order one 
lower than the sun of ohe orders of the constituent row matrices, because 
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in combining these constituent rov matrices the first element of sub- 
matrix [j^oJ added to the last element of . • These sub- 

matrices are, in turn, defined by 




(103) 




and the row matrix [^^oj defined by 




(l)‘ 


r C / CC: 


■hil 


cc- 


f I cC 


cCi 




|l 0 0 . . . (^ (105) 




The quantities 


(Sfel 


are the same as those used in the derivation 


0 


of the matrix in equation (70) 


Division of equation (lOl) by yields 


■'°o 


1 ^ II, 

2Ct= I 
-■So 




(s) 


(106) 


and subseouent premultiplication of this equation by the column matrix 

r_cci (5) 


^r'^L, 


'°0 


w 


yields 
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t 



( 107 ) 


This expression may now he used to eliminate 
tion (96) with the result that 


CC7 

J — 




( 5 ) 


from ecua- 



(108) 


where the aileron reversal matrix 



is defined as 



(109) 


The value of q at aileron reversal is the lowest, real positive value 
of qj^ which satisfies eq-uation (108). 
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TA2LE 1 

GEOHSTTFIC, SIr«UC7URAL, AKB A£EOIT:'NAMIC PASAMETZERS 
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TA3XJE; 2 - Concluded 
SijMMA3Y OF BESULSS 


(t) Supersonic flow 


Wing 

nunber 

Aileron 
conf igu- 
raulon^ 

°T3 


q 

i 

y* 

a 

Po 



■ 

(3) 

-h.56 

-I17.5 

5.55 

0 

3.0 

6.0 

1.000 

.823 

.747 

0.444 

.394 

.363 

0.450 

.371 

.326 

l.COO 

.648 

.494 

1.000 

.250 

-.029 

0.1397 

.0538 

-.0081 

2 

(3) 

-5.75 

Complex 

5.65 

0 

3.0 

6.0 

1.000 

.902 

.870 

.hill 

.412 

.392 

.450 

.436 

.430 

1.000 

-751 

.634 

1.000 

-293 

-.027 

.1397 

.0545 

-.0060 

3 

(3) 

-12.70 

15.8 

5.74 

0 

3.0 

6.0 

1.000 

1.102 

1-303 

.444 

.458 

.472 

.450 

.402 

.357 

1.000 

1.105 

1.326 

1.000 

.542 

-.049 

-1397 

.0686 

-.0052 

it 

(3) 

6.03 


10.85 

0 

3.0 

1.000 

1.427 

.444 

.523 

.450 

.337 

l.OCO 

1.902 

1.000 

1.584 

.1397 

.1163 

5 

(3) 

7.61 



0 

3.0 

1.000 

1-193 

.444 

.493 

.456 

.442 

1.000 

1.494 

1.000 

1.288 

.1397 

.1204 

6 

(3) 

-15.4^ 

18.9 

8.74 

0 

3-0 

6.0 

1.000 

.924 

.964 

.444 

.432 

.432 

.450 

.409 

.391 

1.000 

.926 

.905 

l.OCO 

.650 

.32S 

.1397 

.0981 

.0506 

■ 

(3) 

-134.05 


8.4o 

0 

3.0 

6.0 

1.000 

.997 

.9S0 

.444 

.437 

.435 

.450 

.440 

.437 

1.000 

.981 

.962 

1.000 

.633 

.277 

.1397 

.0901 

.0402 

8 

(3) 

-8.29 

-190.1 

7.10 

0 

3.0 

6.0 

1.000 

.861 

.768 

.444 

.402 

.376 

.450 

.343 

.277 

1.000 
■ 739 
.588 

1.000 

.410 

.084 

.1397 

.0775 

.0200 

9 

(3) 

8.01^ 

n 

m 



.444 

.501 

.450 

-310 

1.000 

1.584 

l.COO 

.983 

-1397 

.0867 

10 

(3) 

4.70 


8.4l 

0 

3-0 

1.000 

1.765 

.444 

.564 

.450 

-125 

i.OOO 

2.615 

1.000 

2.001 

•1397 

.3069 

B 

(1) 

4.70 


6.60 

0 

3-0 

1.000 

1.765 

.444 

.564 

.450 

.125' 

1.000 

2.615 

1.000 

1.633 

.2200 

.1374 

11 

(3) 

GS 


8.50 

0 

3.0 

6.0 

1.000 

1.000 

1.000 

.444 

.444 

,444 

.450 

.450 

.450 

1.000 

1.000 

1.000 

1.000 

.647 

.294 

.1397 

.0612 

.0278 

11 

(1) 

CO 


9.60 

0 

3-0 

1.000 

1.000 

.444 

.444 

.450 

.450 

1.000 

1.000 

1.000 

.688 

.2200 

.1025 

12 

(3) 

-4.22 

-39.3 

■ 

0 

3-0 

6.0 

L.OCO 

.620 

.732 

.444 

.382 

-350 

.450 

.275 

.185 

1.000 

.623 

.465 

1.000 

.215 

-.057 

■ 1397 
.0482 
-.0172 

12 

(1) 

-4.22 

-39.3 

8.57 

0 

3.0 

1.000 

.820 

.444 

.382 

.450 

.275 

i.OOO 

.623 

1.000 

.355 

.2200 

.1255 


^Aileron configurations: 

(1) 50— percent-senispan outboard aileron, 20 percent chord. 

(2) 30-percent-senispan ounboard. aileron, 100 percent chord 

(3) 30-?e rcen w-se3ispan outboard aileron^ 20 percent chord- 
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TABLE 3 

THS INT3GRATi:?G l^ATEICSS 


(a) y*B = 0.3 


Ml 


B 

.1 

.2 

• 3 

B 

0 . 16667 

0.66667 

0.16667 

E 

-.c 4 i 67 

•33333 

.20833 

.3 

0 

0 

0 




y* 

.30 

.hk 

.58 

.72 

.86 

.30 

0.06667 

0.25000 

0.18333 

0.17895 

0.27000 

.kk 

-.01667 


.20000 

. 17895 

.27000 

.58 

0 

-.01667 

. 11667 

.17895 

.27000 

.72 

0 

0 

- .01667 

.09562 

.27000 

.86 

0 

0 

0 

-.03771 

. 18667 


[ 13 ] 



i-l 

.2 

• ^ 

B 

0 

0.33333 

0 . 16667 

E 

-.01042 

.06250 

.07292 


0 

0 

0 


C1I-] 

o 


y* 

.30 

.Lh 

.58 

.72 

.86 

• 30 

0 

0 

• 

0 

§ 

0 

0.07833 

0 . 10052 

0.22776 

.44 

- 

00167 

.00500 

.04000 

.06473 

.17376 

.58 

0 


.00500 

.02894 

.11976 

.72 

0 

0 

-.00167 

-.00185 

.06576 

.86 

0 

0 

0 

-.00431 

.01676 


B 

.1 

.2 

. 3 

B 

0 

0 

0 


.20833 

•33333 

-. 04 i 67 

■ 

. 16667 

. 6666"^ 

. 16667 


&]"o 


y* 

.30 

.44 

.58 

.72 

.86 

• 30 

0 

0 

0 

0 

0 

.44 

.08333 

.13333 

-.01667 

0 

0 

.58 

.03333 

.21667 


-.01667 

0 

.72 

.08333 

.21667 

.20000 

. 11667 

-.01667 

.86 

.08333 

. 21667 

.20000 


. 11667 
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TA3EE 3 - Continued I 
THE INTEGRATING MATRICES 


(a) y*^ = 0.5 



y* 

.1 

.2333 

.3667 

■ 5 

.1 

O-lllll 

0 . 41667 

0.33333 

0 . 13889 

.2333 

-.02778 

. 1944 Ji 

. 36111 

.13889 

.3667 

0 

-.02778 

.22222 

.13889 

.5 

0 

0 

0 

0 


&’]o 


.y* 

.5 

.625 

.75 

.875 

.5 

0.08333 

0.31250 

0.20286 

0 . 33750 


-.02083 

.14583 

.22369 

.33750 

.75 

0 

- .02083 

.11953 

.33750 

HS 

0 

0 

-.04714 

.23333 


Ml 


y * 

. 1 

.2333 

.3667 

.5 

, 1 

0 

0 . 12500 

0.25000 

0 . 12500 

.2333 

-.00^^63 

.01389 


.07870 

.3667 

0 

-.00463 

.02778 

.03241 

.5 

0 

0 

0 

0 




y* 

.5 

.62p 

.75 

.875 

.5 

0 

0.07031 

0.09854 

0.27150 

.625 

-.00260 

.00781 

.04522 

.18713 

.75 

0 

-.00260 

-.00289 

.10275 

.875 

0 

0 

- .00673 

.02619 




.1 

.2333 

.3667 

.5 

.1 

0 

0 

0 

0 

.2333 

.15889 

.22222 

-.02778 

0 

.3667 

.13889 

. 36111 

.19444 

-.02778 

m 

.13889 

.33333 

. 4 i 667 

.11111 


y* 

.5 

.625 

.75 

.875 

.5 

0 

0 

0 

0 

.625 

.10417 

.16667 

-.02083 

0 

.75 

. 10 417 

.27083 

.14583 

-.02083 

.875 

. 10417 

.27083 

.25000 

.14583 
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1 TA3LS 3 - Continued 
THE INTEGRATING MATRICES 


(c) y*^ = 0.7 
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TABLE 3 - Concluded 
TEZ INTEGRATING MATRICES 

(<3.) = 1 or 0 (ordinary vlngs) 


y* 

.10 

■ 25 

.■to 

-55 

.70 

.85 

.10 

0. 05550 

0.20833 

0.15278 

0.16667 

0.14913 

0.22500 

.25 

-.01389 

.09722 

.16667 

.16667 

. 14913 

.22500 

.1(0 

0 

-.01389 

.09722 

.16667 

.14913 

.22500 

.55 

0 

0 

-.01389 

.09722 

.14913 

.22500 

.70 

cr 

0 

0 

-.01389 

.07968 

.22500 

.85 

0 

0 

0 

0 

-.03143 



-.00116 


.25 


0.03125 


. 003^^7 


-.00116 


0 


0 


0 


Li: 


• to 


0.05i^^^0 


.02778 


.003^7 


- .00116 


.55 

.70 

.85 

0.08333 

0 .09466 

0.19567 

.05556 

.06980 

. 15817 

.02778 

.04495 

.12067 

.00347 

.02010 

.08317 

-.00116 

-.00129 

.04567 

0 

- .00299 

.01164 


.10 

.25 

.to 

0 

0 

0 

.06914 

.11111 

-.01389 

.06944 


.09722 

.06944 

.18056 

.16667 

.06944 

. 18056- 

.16667 

.06914 

.18056 

. 16667 


-.01389 


.09722 


.16667 


.16667 


.70 

.85 

0 

0 


0 


-.01389 


.09722 


.16667 


0 


-.01389 


.09722 
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(a) M- wings. 

Figure 2.- The spanwisc variations of stiffness. 


■•ji 





Dimens. 
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(d) Swept wings and unswept wing. 
Figure 2.- Concluded. 
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eciion momenf arms, e/ and 62. 


.2 .3 .4 .5 .6 .7 .8 

Dimensionless spanwise ordinate, y* 


(■fa) W-wings. 


Figure 3 .- Continued. 
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( c ) A-wings and inverted A-wings . 
Figure 3.- Continued. 




Secii 



Dimensionless span wise ordmaie, 


(d) Swept wings and unswept wing. 
Figure 3.- Concluded. 
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(a) Lift distributions due to unit airplane angle of 

attack (subsonic). 

(b) Lift distributions due to iinit effective aileron 

deflection (subsonic). 


Figure 4 .- Spanwlse lift distributions for wing 1. 


^M-wingj y*g = 0.3.) 
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Dimensionless spaznvlse ordinate , y* 

(c) Lift distri’outions due to unit airplane angle of 

attack ( supersonic ) . 

(d) Lift distributions due to unit effective aileron 

deflection ( supersonic ) . 


Figure il-.- Concluded 
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(a) Lift distributions due to \anit airplane angle of 

attack ( subsonic ) . 

(b) Lift distributions due to unit effective aileron 

deflection ( subsonic ) . 

Figure 5-- Spanwise lift distributions for -wing 2. ^M-wing; y*g = 0.5. 


Spanwise loadlin; coefticlent, — Spajiwjse loading coefficient 
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Dimensionless spanwise ordinate, y* 


(c) Lift distributions due to unit airplane angle of 

attack ( supersonic ) . 

(d) Lift distributions due to unit effective aileron 

deflection ( supersonic.) . 


Figure 5 .- Concluded 
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(a) L,ift distributions due to -unit airplane angle of 

attack ( subsonic) . 

(b) Lift distributions due to unit effective aileron 

deflection (subsonic). 

Figure 6.- Spanwise lift distributions for wing 3. ^M-wingj y*g = 0.7.) 
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(c) Lift distributions due to unit airplane angle of 

attack ( supersonic ) . 

(d) Lift distributions due to unit effective aileron 

deflection ( supersonic ) . 


Figure 6.- Concluded 
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(a) Lift distributions due to unit airplane angle of 

attaclc ( subsonic ) . 

(b) Lift distributions due to unit effective aileron 

* deflection (subsonic). 

Figure 7 .- Spanwise lift distributions for wing 4 . (w-wing; = 0.3.) 
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(c) Lift distrituticr-B due tc unit airplane angle of 

attack ( supersonic ) . 

(d) Lift distributions due to unit effective aileron 

deflection ( supersonic) . 


Figure 7 .- Concluded 
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Dimensionless spair<vlse ordinate , y* 


(a) Lift distributions due to unit airplane angle of 

attack ( subsonic ) . 

(b) Lift distributions due to unit effective aileron 

deflection ( subsonic ) . 


Figure 8.- Spanwise lift distributions for wing 5. ^W-wing; = 0.5* ) 
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(c) Lift distributions due to unit airplane angle of 

attack (supersonic). 

(d) Lift distributions due to unit effective aileron 

deflection (supersonic). 


Figure 8.- Concluded 


Spanwlse loading coefficient, i- Spaiwlse loading coefficient 
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jDlir.ensionless spacwlse ordinate, y* 


(a) iiift distributions due to unit airplane angle of 

attack (subsonic). 

(b) Lift distribu'cions due to unit effective aileron 

deflection (subsonic). 

Figure 9.- Spanwise lift distributions for wing 6. ^W-wing; y*-g = 0.7.) 
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(c) Lift distributions due to unit airplane angle of 

attack ( supersonic ) . 

(d) Lift distributions due to unit effective aileron 

deflection ( supersonic ) , 


Figure 9-- Concluded 
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Spanwise loading coefiicient, - — — Spanwise loading coeMcient 
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(c) Lift distritutlons due to unit airplane angle of 

attack ( supersonic ) . 

(d) Lift distributions due to unit effective aileron 

deflection ( supersonic ) . 


Figure 10.- Concluded. 


nn 




Spanwlse loading coefficient, Spanwlse loading ooefflclent 
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0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

Dimensionless spanwise ordinate, 


(a) Lift distributions due to unit airplane angle of 

attack (subsonic). 

(b) Lift distributions due to unit effective aileron 

deflection (subsonic). 

Figure 11.- Spanwlse lift distributions for wing 8. ^A-wlng; y*g = 


0.7.) 
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(c) Lift distributions due to unit airplane angle of 

attack (supersonic). 

(d) Lift distributions due to unit effective aileron 

deflection ( supersonic ) . 


Figure 11.- Concluded 
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(a) Lift distributions due to unit airplane angle of 

attack (subsonic). 

(b) Lift distributions due to ■unit effective aileron 

deflection ( subsonic ) . 


Figure 12.- Spanwise lift distributions for wing 9 - 

y*B = 0 . 7 .) 


(inverted A-wlng; 



Spanwlse loading ooeHlolent, ,=^ Spaiwise loading coefficient, 
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.8 .9 


Dimensionless spanwlse ordinate, y* 


(c) Lift dlstrllmtions due to unit airplane angle of 

attack ( supersonic ) . 

(d) Lift distributions due to unit effective aileron 

deflection ( supersonic ) . 


Figure 12.- Concluded. 
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(a) Lift distributions due to iinit airplane angle of 

attack ( subsonic ) . 

(b) Lift distributions due to unit effective aileron 

deflection ( subsonic ) . 

Figure 13.- Spanwise lift distributions for wing 10. ( Sweptf orward wing.) 
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(c) Lift distributions due to unit airplane angle of 

attack ( supersonic ) . 

(d) Lift distributions due to unit effective aileron 

deflection (supersonic). 

Figure 13.- Concluded. 
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(a) Lift distributions due to lanit airplane angle of 

attack ( subsonic ) . 

(b) Lift distributions due to \mit effective aileron 

deflection ( subsonic ) . 

Figure l4.- Spanwlse lift distributions for wing 11. (Unswept wing.) 


Spanwise loading coefficient, — Spanwise loading coefficient 
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(c) Lift distributions due to unit airplane angle of 

attack ( supersonic ) . 

(d) Lift distributions due to imit effective aileron 

deflection ( supersonic ) . 


Figure l4.- Concluded. 



Spanwise loading codfflclent, Spanwiae loading coefficient 
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(a) Lift distributions due to unit airplane angle of 

attack (subsonic). 

(b) Lift distributions due to \mit effective aileron 

deflection (subsonic). 


Figure 15.- Spanwise lift distributions for wing 12. (Sweptback wing.) 
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(c) Lift distributions due to \mit airplane angle of 

attack ( supersonic ) . 

(d) Lift distributions due to unit effective aileron 

deflection ( supersonic ) . 


Figure 15.- Concluded 
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Dynamic-pressure poramefer, q 


Figure 17.- Aeroelastlc characteristics of wing 2. ^M-wing; y*g = 0.5.) 
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Dynamic-pressure parameter, q 

Figiire l8.- Aeroelastlc characteristics of wing 6. ^W-wlngj y*g = 0 . 7 .) 
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Figxire 20.- Aeroelastlc characteristics of wing 8. f A 
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Subsonic 




1 . 2 ^ 



Dynamic-pressure parameter, q 


Figure 21.- Aeroelastic characteristics of wing 11, (Unswept wing.) 
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(a) M-wings. 

Figure 23.- 'Variation of aeroelastic characteristics vith ving "break 

position; q = 3-0. 
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Spanw/se ordmate at v/mg break, 


(b) W-wings. 
Figure 23.- Continued 
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Figure 2h. - Concluded. 
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(a) Swepfbac2c inner pa 

(b) Swept forward inner p 

Figure 25.- Combinations of and y*j^ 

zero aerodynamic-center shift 






Angle of sweep of outer wing panel, yi 
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